Optical transitions in highly-charged californium ions with high sensitivity to 

variation of the fine-structure constant 
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We study electronic transitions in highly-charged Cf ions that are within the frequency range of 
optical lasers and have very high sensitivity to potential variations in the fine-structure constant, 
q. The transitions are in the optical despite the large ionisation energies because they lie on the 
level-crossing of the 5/ and 6p valence orbitals in the thallium isoelectronic sequence. Cf 16+ is a 
particularly rich ion, having several narrow lines with properties that minimize certain systematic 
effects. Cf 16+ has very large nuclear charge and large ionisation energy, resulting in the largest 
Q-sensitivity seen in atomic systems. The lines include positive and negative shifters. 

PACS numbers: 06. 30. Ft, 31. 15. am, 32.30.Jc 



I. INTRODUCTION 

In this work we present calculations of transitions in 
highly-charged californium that could form the reference 
for an optical atomic clock with very strong sensitivity 
to variation of the fine-structure constant, a — e 2 /Tic. 
Our work is motivated by recent astronomical studies of 
quasar absorption spectra that indicate a spatial gradient 
in values of a across cosmological distances [l|, Q . The re- 
sults were taken using around 300 spectra covering most 
of the sky, observed at two telescopes: the Very Large 
Telescope in Chile Q and Keck Telescope in Hawaii pj- 
5]. The telescopes independently agree on the direction 
and magnitude of the gradient (dipole) , which is signifi- 
cant at 4.2a for the combined sample of both telescopes. 

The cosmological dipole in a might be confirmed by 
terrestrial studies, since the solar system is moving with 
respect to the cosmic microwave background (the pre- 
sumed frame for the ce-dipole), and therefore should be 
moving from a region of the Universe with smaller values 
of a to one with larger values @. In particular, the ex- 
pected rate of change in a today would be of order a/ a ~ 
10 _18 yr _1 . This is significantly smaller than the best 
current terrestrial limits, a/a = (— 1.6±2.3)xl0 _17 yr _1 , 
which comes from comparison of Al + and Hg+ atomic 
clocks 0- If measured at the same level of accuracy, 
the transitions proposed in this work would allow an im- 
provement on this limit by a factor of 23. Because the 
transitions have narrow natural line widths and reduced 
systematics, the improvement could be even larger. 

We parametrize the sensitivity of an atomic transition 
to potential variation in a by the quantity q defined by 



dx 



(1) 



where x — a 2 / a\ — 1 is the fractional change in a 2 from 
its current value a\ , and q and u are measured in atomic 
units of energy. In the A1+ and Hg+ comparison, the 
A1+ clock is an "anchor" (relatively insensitive to co- 
variation) while the mercury clock has a strong sensitivity 
of q = —52 200 cm -1 [8]. An approximate formula for the 



q value of a single energy level (E n = —I n where /„ is 
the ionisation energy of the level) with effective principal 
quantum number v and angular momentum j is ficil ] 



Qn ~ In 



(Za) 2 



v{j + 1/2) ' 



(2) 



where Z is the nuclear charge. The transition will have 
a sensitivity to a- variation that is the difference between 
the q values of the levels involved. Therefore the best 
transitions will maximise the difference of v and j be- 
tween the levels and will come from heavy ions. 

Equation [2] shows that transitions in highly charged 
ions (HCIs) can have much larger q values since they have 
much larger ionisation energies. Unfortunately they gen- 
erally also have much larger transition energies, putting 
them outside the range of optical lasers and making them 
unsuitable for use in high-precision clocks. However, due 
to configuration crossing, some HCIs can have optical 
transitions between levels with different principal quan- 
tum numbers, and these could become reference transi- 
tions for optical clocks with the highest q values seen in 
atomic systems [ltj [H| . 

In this paper we present calculations for the 5/ - 6pi/2 
crossing that occurs in the thallium and lead isoelectronic 
sequences (with one valence electron and two valence 
electrons, respectively). The crossing occurs at Z = 98 
- californium — which also happens to be one of the 
last relatively stable ions in the periodic table, having 
isotopes with half-lives of several hundred years. The 
five isotopes shown in Table U allow for choice in exper- 



TABLE I. Half-lives and symmetries of long-lived californium 
isotopes [121 ]. 



Isotope 


r 


Half-life (yr) 


248 Cf 


0+ 


0.914 


249 Cf 


9/2" 


351 


250 Cf 


0+ 


13.08 


251 Cf 


1/2+ 


900 


252 Cf 


0+ 


2.645 



2 



FIG. 1. Dirac-Fock energies of the 6pi/2 (diamonds, dashed 
line), 6P3/2 (crosses, dot-dashed line), and 5/ (circles, solid 
line) levels in the thallium isoelectronic sequence with increas- 
ing nuclear charge. The inset shows an enlarged view of the 
crossing region. 
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iments: isotopes with an odd number of nucleons have 
a greater half-life, but isotopes with an even number of 
nucleons exclude the possibility of hyperfine structure. 
We find that Cf 16+ has several transitions that may be 
suitable for studies of a-variation, including those with 
the largest q values yet found. 



II. METHOD 

To find the 6p - 5/ crossing, we start with neutral thal- 
lium. In thallium, the 5/ orbital energies lie above the 6jj 
orbitals, whereas in the large Z limit the 5/ levels should 
be more tightly bound than 6p levels (E$f « E^p for 
hydrogen-like ions). Therefore we expect a level crossing 
at some Z > 81, where an ion may have optical transi- 
tions between these two orbitals. Fig. Q] shows the Dirac- 
Fock energies of the 6pi/ 2 , 6j?3/2, 5/5/2 and 5/7/2 orbitals 
as a function of nuclear charge Z. Due to the large fine- 
structure splitting of the Qpi/2 and 6p 3 / 2 subshells, there 
are two possible crossing points we can explore here, one 
for 6p 3 / 2 near Z — 93 and one for 6p!/ 2 near Z = 98. The 
crossing near Z — 98 is more attractive for studying co- 
variation for two reasons. Firstly, the nuclear charge and 
ionization energy are larger. Second, since a-sensitivity 
is due to relativistic effects that occur near the origin, 
the 6p!/ 2 orbital has larger q than the 6p 3 / 2 orbital since 
the former has a lower Dirac-spinor component of S1/2 
symmetry, which is large near the origin. This is seen in 
@ by the factor + 1/2): the difference in q due to 
this factor is larger for P1/2 and fs/2- Because the 6pi/2 
level is highly sensitive to a-variation while the 5/ levels 
are not, we expect a large q value for a transition between 
these levels. 

We have performed full-scale ab initio calculations for 
Cf 1,+ (at the crossing point of Fig. [1} and the two- 



TABLE II. Low-lying levels of Cf 17+ (all have odd parity). 
Energy calculations are presented relative to the 5/5/2 ground- 
state using only Dirac-Fock (DF) and including MBPT (DF 
+ E). The q values were calculated using DF + E. 



Configuration 


J 


Energy ( 


cm x ) 


q (cm ) 






DF 


DF + E 




5/ 


5/2 











6p 


1/2 


8447 


17889 


-449750 


5/ 


7/2 


20447 


21755 


17900 


6p 


3/2 


233514 


241970 


-115650 



valence-electron equivalent, Cf 16+ . We use the combined 
configuration interaction and many-body perturbation 
theory method (CI + MBPT), presented in detail in [l3[ 
(see also [HI). We begin with Dirac-Fock for closed shells 
of Hg; this corresponds to V^ 1 for the single-valence- 
electron case, Cf 17+ , and V* -2 for the Cf 16+ ion. From 
the frozen-core potential we generate a set of around 
40 B-splines in each wave up to / = 6. These form a 
"complete" set of virtual orbitals with which we calcu- 
late MBPT corrections, E, to second-order in the resid- 
ual Coulomb interaction. For Cf 16+ we perform a CI 
calculation including all two-electron excitations to the 
virtual orbitals IQspdf. The addition of g-wave orbitals 
to the CI were found to make little difference to energy 
levels and q values. The q values were obtained by vary- 
ing x (Eq. [T]) in steps of 0.01 and taking the gradient of 
transition frequency with respect to x. 

We have also calculated some important transi- 
tion rates (reduced matrix elements and Einstein A- 
coefficients) using a relativistic formalism (see, e.g. fl5|). 
Random phase approximation corrections to the matrix 
elements were not included in this work since the uncer- 
tainty in the rates is dominated by uncertainty in the 
transition energies, which have not been measured. 



III. RESULTS AND DISCUSSION 

Our calculated energy levels and q values for Cf 17+ 
are presented in Table [TTJ For the purposes of measur- 
ing a-variation using atomic clocks, the most interesting 
transition is from the 5/ 5 / 2 ground state to 6pi/ 2 , with an 
energy interval of cj = 17889 cm -1 , q for this transition 
is over 8 times larger than the Hg + clock transition used 

in a. 

In Cf 16+ , presented in Table IIII1 there are more states 
in optical range, arising from the greater number of an- 
gular momentum combinations available. Note that the 
levels marked A, A* in Table [TTTI are heavily mixed in the 
CI calculation, resulting in a dominant contribution from 
5/ 2 (63% and 51% respectively) in both, while in the CI 
+ E calculation the first state is 96% 5/ 2 and the second 
state is 63% 5/ 6p. For simplicity we have simply labelled 
the level with the largest 5/ 6p contribution as the 5/ 6p 
state. 
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TABLE III. Calculated energy levels, g-factors, and q values 
for low-lying levels of Cf 16+ (all have even parity) relative to 
the ground state 5/6p (J = 3). Energies are calculated using 
only configuration interaction (CI) and including many-body 
perturbation theory effects (CI + E). The ID column just 
provides a convenient label reference in the text. 



Config. J ID g Energy (cm 1 ) q (cm ) 











CI 


CI + E 




5f6p 


3 


Gl 


0.8299 











6p 2 





G2 




-7429 


5267 


-370928 


5/6p 


2 




0.8482 


7313 


6104 


106124 


5/ 2 


4 


A 


0.8535 


28746 


9711 


414876 


5/6p 


4 


A* 


1.0481 


21415 


24481 


162126 


5/ 2 


2 




0.7532 


38674 


24483 


354444 


5f6p 


3 




1.1776 


23979 


25025 


59395 


5/ 2 


5 




1.0333 


43097 


29588 


451455 


5/ 2 


3 




1.0954 


50953 


37467 


393755 


5/ 2 


4 




1.1197 


53229 


42122 


319216 


5/ 2 


6 




1.1371 


57220 


44107 


459347 


5/ 2 





B 




68192 


51425 


380986 


5/ 2 


2 




1.1672 


67267 


51471 


446045 


5/ 2 


4 




1.1198 


69475 


58035 


461543 


5/ 2 


1 


C 


1.5000 


75018 


58132 


449977 


5/ 2 


6 




1.0296 


78739 


63175 


460416 


5/ 2 


2 




1.2672 


89580 


75041 


465293 


5/ 2 









127521 


114986 


446376 


5/6p 


3 




0.9765 


211414 


212632 


323435 


6p 2 


1 




1.4963 


198879 


213864 


-113277 



Our Cl-only calculations showed the 6p 2 (J = 0) level 
to be the ground state, but adding MBPT corrections 
changes the level ordering such that 5/ Qp ( J = 3) is the 
ground state. Actually, Cf 16+ can be considered to have 
two ground states, since the decay from the metastable 
6p 2 (J = 0) (G2 in Table!]!]) to the ground state (Gl) has 
a lifetime greater than that of the nucleus itself. Table ITVl 
lists calculated matrix elements and strengths for some 
transitions of interest. 

The two electron transitions between the G2 
metastable state and the 5/ 2 states give maximal val- 
ues of q: up to around ~ 830 000 cm -1 . Among these is 
the transition to 5/ 2 (J = 1) (C in Table IHlT) with en- 
ergy cj = 58 132 — 5267 = 52 865 cm -1 , which has a very 
high branching-ratio back to the G2 "ground" state. This 
level therefore potentially provides a method to "recycle" 
from Gl back to G2, although it should be noted that 
the Gl — > C transition is rather weak. 

Another very interesting potential reference transi- 
tion is the G2 (J = 0) — > B (J = 0) transition at 
uj = 46158 cm -1 , which is strongly forbidden but could 
be opened using Stark amplitude or hyperfine mixing of 
state B (J = 0) with C (J = 1). Such a transition 
would be very narrow and have strongly reduced system- 
atic shifts, e.g. electric quadrupole, AC Stark, Zeeman 
shifts. It may, however, be too weak to excite by usual 
optical lasers. 

All of the transitions discussed so far are positive 



shifters: the transition frequency increases with increas- 
ing a. It is also possible to find negative shifters in Cf 16+ , 
for example the transition between Gl and G2 is a strong 
negative shifter (assuming that the ordering of levels has 
been calculated correctly). However, this transition is 
extremely weak, and in practise may only occur via level 
mixing using a strong laser. A negative shifter which may 
be more useful is from the 5/ 2 (J = 4) metastable state 
(A in Table |nTl lifetime ~ 10 -1 s) via Ml transition to 
one of the 5/ 6p states above it. The larger of these has 
q = -355 000 cm -1 . 

IV. SYSTEMATICS AND OPPORTUNITIES 

HCIs have some interesting features that are worth 
mentioning here. Firstly, electric dipole matrix elements 
are much smaller for HCIs than in neutral atoms since 
the El matrix element ~ (r) ~ (ao/Z c g) where ao is the 
Bohr radius and Z e s « Z lon + 1 is the effective nuclear 
charge: the charge that the valence electron sees. Since 
the spacing between El levels in HCIs is larger by a factor 
~ Z 2 ff , the static polarizability — and hence black-body 
radiation shift — of HCIs is reduced compared to near- 
neutral ions by a factor ~ 1/Z^ S . 

The hyperfme structure in heavy HCIs is much larger 
than in neutral atoms, scaling as u>hfs ~ ZZ 2 S . The 
rate of Ml transitions within each hyperfine multiplet 
will scale as w^y s , which means that the lowest hyperfine 
state will be produced in reasonable time (order of a sec- 
ond). In californium, the hyperfine splitting of an s-wave 
or Pi/2-wave valence electron will be very sensitive to co- 
variation because of the large Z. We define the fractional 
(relative) sensitivity K IC \ by Suihfs/uJhfs = K re i5a/a. 
Using formulas presented in Refs. [1(1 G3] we obtain 
Kiel — 5.33. Therefore, the hyperfine transitions form 
another positive shifting transition that can be used to 
place limits on a- variation. 

V. CONCLUSION 

In this paper, we used the CI + MBPT method to 
calculate the energy levels for Cf 16+ and Cf 17+ highly 
charged ions. These ions were chosen because they lie 
at the 5/ - 6pi/2 crossing point on the thallium isoelec- 
tronic sequence, which allows for optical transitions be- 
tween different configurations from the ground state. Our 
calculations have identified several transitions in Cf 16+ 
that have the largest q values ever seen in such atomic 
systems, and include several positive shifters (with q up 
to ~ 830 000 cm -1 ) and negative shifters (e.g. 5/ 2 (J = 
4) -> 5/6p (J = 3) with q = -355 000cm -1 ). A com- 
parison of clocks using these reference transitions would 
have a total sensitivity Aq = q + — q_ w 1 200 000 cm -1 , 
around 23 times more sensitive than the Hg + clock and 
Al + clock comparison used to obtain the best current 
laboratory limit on cv-variation. 
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TABLE IV. q values, squared reduced transition matrix elements, S, and corresponding Einstein ^-coefficients for transitions 
between selected states i and / in Cf 16+ . The included 6p 2 has longer lifetime than the nucleus. 



Config. 



cm- 1 ) 


Config. 


Jf 

J 


E f (cm -1 ) 


Pit (cm 1 ) 


S 


at. An (s x ) 










5267 


-370928 


0.92401m3 


4.3519xl0~ ls 




5f 6» 


2 


6104 


106124 


0.14553mi 


0.89281 












0.35612 B2 


0.00033798 




5f 2 


4 


9711 


414876 


0.16895mi 


4.1731 










11 


1.6768b2 


0.016216 




5/6p 


4 


24481 


162126 


2.5836mi 


1022.5 












0.20237e2 


0.1993 




5f 2 


2 


24483 


354444 


0.0041938mi 


1.6601 












0.079461 B2 


0.078285 




5/6p 


3 


25025 


59395 


0.071521 M i 


30.235 










n 


0.0022665 B 2 


0.0024914 




5/ 2 


5 


29588 


451455 


0.0011599 B2 


0.0029459 




5/ 2 


3 


37467 


393755 


0.063566 M i 


90.182 










11 


0.016359 B2 


0.13528 




5/ 2 


4 


42122 


319216 


1.1631 M i 


2344.7 










n 


0.25404 B 2 


3.7727 




5/ 2 


2 


51471 


446045 


0.0082138mi 


30.211 










11 


0.016684 B2 


0.6750 




5/ 2 


4 


58035 


461543 


0.13984 M i 


737.27 










n 


0.028702 B2 


2.1162 




5/ 2 


1 


58132 


449977 


0.0023199 B2 


0.17249 




5/ 2 


2 


75041 


465293 


0.003621 5 M i 


41.278 










11 


0.0022201i52 


0.59165 


5267 


5/6p 


2 


6104 


477052 


0.54123i52 


2.4845xl0~ 8 




5/ 2 


4 


9711 


785804 


0.24444 S4 


8.5359xl0" 24 




5/ 2 


2 


24483 


725372 


0.019556 B2 


0.0057377 




5/ 2 


2 


51471 


816973 


0.020833s2 


0.49127 




5/ 2 


1 


58132 


820905 


0.1195mi 


476.2 




5/ 2 


2 


75041 


836221 


0.00071052 B2 


0.13159 


6104 


5/ 2 


4 


9711 


308752 


0.1839 E2 


1.2569xl0~ 5 




5/6p 


4 


24481 


56002 


0.028394 B2 


0.0066653 




5/6p 


3 


25025 


-46729 


3.1428mi 


574.24 




5/6p 


3 


25025 


n 


0.033026 B2 


0.0089701 




5/ 2 


1 


58132 


343853 


0.0082881mi 


31.485 




5/ 2 


1 


58132 




0.043495 B2 


1.8571 


9711 


5f 6p 


4 


24481 


-252750 


1.605mi 


139.5 










11 


0.052382 B 2 


0.0041243 




5/ 2 


2 


24483 


-60432 


1.5203 B2 


0.11978 




5/6p 


3 


25025 


-355481 


0.096413mi 


9.3407 












0.0016877 B2 


0.00015922 




5/ 2 


5 


29588 


36579 


8.3595 mi 


1771.0 








11 


0.023016 B2 


0.0079994 






5/ 2 


3 


37467 


-21121 


0.06325mi 


36.484 












0.055057b2 


0.10159 




5/ 2 


4 


42122 


-95660 


0.39244 M i 


360.42 












0.26664i52 


1.0681 




5/ 2 


6 


44107 


44471 


0.00026242 B2 


0.0014151 




5/ 2 


2 


51471 


31169 


0.035908b2 


0.51073 




5/ 2 


4 


58035 


46667 


0.0011045mi 


3.362 










11 


0.0033787i52 


0.099715 




5/ 2 


6 


63175 


45540 


5.15X10" 6 B2 


0.00025212 




5/ 2 


2 


75041 


50417 


0.015122 B 2 


2.0155 


24481 


5/6p 


3 


25025 


-102731 


0.99502mi 


0.0043219 


24483 


5/ 2 


1 


58132 


95533 


0.024099mi 


24.766 










11 


0.18074 B 2 


0.87324 


25025 


5/ 2 


1 


58132 


390582 


0.084098b2 


0.37462 


37467 


5/ 2 


1 


58132 


56222 


1.1076 E 2 


0.46748 


51425 


5/ 2 


1 


58132 


68991 


1.6793mi 


13.671 


51471 


5/ 2 


1 


58132 


3932 


1.0698mi 


8.5301 










11 


0.58299E2 


0.00085652 



5/6p 



6p 2 



5/6p 



5/ 2 



5/6p 
5/ 2 



5/6p 
5/ 2 
5/ 2 
5/ 2 
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Trapping and cooling of HCIs remains a difficult exper- 
iment, however electron-beam ion trap technology con- 
tinues to improve [lsI - EII ] , and we hope that the potential 
benefits of HC1 clocks will continue to motivate further 
studies. Using californium certainly adds another layer of 
complexity to the experiment since it doesn't occur natu- 
rally and must be produced, for example, at accelerators 
such as GSI and LBNL. 
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Appendix A: Transition matrix elements 

The CI+MBPT method produces wavefunctions \I) 
which are linear combinations of Slater determinants of 
orbitals (we use atomic units ft = e = m e = 1) 

^ nKm = r \ia Q ' (A1) 

To calculate the relativistic electric (EJ) and magnetic 
(MJ) multipole reduced matrix elements between two 
wavefunctions ipi and tpf we use the following formulae 
(see, e.g. [HJ): 




where 



(n lKl \\q^\\n fKf ) = (2 J ± 1} " (- Ki \ \C J \ \ Kf ) j ^(kr) [fi(r)gf(r)+ gi (r)f f (r)} dr 



(n 4 Ki||gS B) l|n/«/> = ^7^ (^11^ 11*/) / Mkr)[fi(r)ff(r)+ 9i (r)g f (r)} + 



(2J+1 



,!! 



h+i(kr) l-j+f lfi(r)9f(r)+ gi (r)f f (r)} + [Mr)g f (r)- gi (r)f f (r)] 1 dr 



In these equations, 



(«i||C J H«/> = (-l) J ' +1/ V(2j 4 + l)(2 J/ + l) x 



-1/2 l/2 ijM+lf + J) (A6) 



where 



1, for x even 
0, for x odd , 



(A7) 



and jj (kr) is the spherical Bessel function of order J with 
argument kr 



— a;/l r. We also calculate Einstein A- 



I 

coefficients using 

. (A) _ (2J + 2)(2J + l)fc 2 - 7 + 1 



(I\\Q { j X) \\F) 



J[(2J+1)!!] 2 



(2 J F + 1) 



(A4) 
(A5) 



(A8) 



where the capital letters replacing previous lower case 
letters indicate that the quantity should be summed over 
the entire many body system, for example 



Qjm = J2(ljM)if a \ a f 

if 



(A9) 
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